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We study the low energy states of the t-J model with two holes on a 32-site lattice with periodic 
boundary conditions. In contrary to common belief, we find that the state with d x 2_ y 2 symmetry 
is not always the ground state in the realistic parameter range 0.2 < J/t < 0.4. There exist low- 
lying finite- momentum p-states whose energies are lower than the d x 2_ y 2 state when J/t is small 
enough. We compare various properties of these low energy states at J/t — 0.3 where they are 
almost degenerate, and find that those properties associated with the holes (such as the hole-hole 
correlation and the electron momentum distribution function) are very different between the d x 2_ y 2 
and p states, while their spin properties are very similar. Finally, we demonstrate that by adding 
"realistic" terms to the t-J model Hamiltonian, we can easily destroy the d x 2_ y 2 ground state. This 
casts doubt on the robustness of the d x 2_ y 2 state as the ground state in a microscopic model for 
the high temperature superconductors. 

PACS numbers: PACS: 71.27. +a, 71.10.Fd, 75.40.Mg 



I. INTRODUCTION 

The t-J model was proposed as a microscopic model to 
describe the low energy physics of high temperature su- 
perconductors. It mimics the doped CuC>2 planes found 
in high temperature superconductors by a system con- 
sisting of mobile holes moving in a spin background on a 
two-dimensional square latticed One possible mechanism 
for superconductivity is the Bose-Einstein condensation 
of "preformed" hole pairs. To see whether this is possi- 
ble, first we have to understand how the holes interact. It 
is well-known that a mobile hole causes frustration in the 
antifcrromagnetic spin background. Longer range inter- 
action between the holes must be mediated .through the 
spin distortion they produce. Various studiesB seem to in- 
dicate that in the "realistic" parameter range J/t ~ 0.3- 
0.4, the holes form bound pairs in the d x 2_ y 2 channel. 
Theoretically this can be understood as resulting from 
the magnon exchange.interaction which selects the d x 2_ y 2 
channel over others. Q Alternatively, the density matrix 
renormalization group (DMRG) approach has been used 
to work out the spin structure in the vicinity of a hole 
pair in real space O It was found that by pairing up at 
a distance of y/2, the two holes can share the spin frus- 
tration and is therefore energetically favorable. This is 
consistent with the numerical observationErErCl that the 
holes have the largest probability of being at a/2 apart. 
Later it was shown that such pairing of holes at v2 arises 
naturally if the wavefunction has d x 2_ y 2 symmetry!! But 
in order to firmly establish the role of d-w&ve hole pairing 
in the theory of superconductivity, one has to understand 
how strong and robust the pairing is. In fact numerical 
studylj has shown that the binding energy of the d chan- 
nel is negative but small, and that the holes may not be 
tightly bound in real space. Recently the question on 
the robustness of d-wave pairing was addressed by using 
the anisotropic t-J z model as the starting point.EI It con- 
cluded that different mechanisms may select states with 



different symmetries and their competition may destroy 
the d-w&ve ground state. This raises doubt on the ro- 
bustness of the d-wave ground state because it is known 
that the t-J model alone is not enough to explain fully 
the hole dynamics of StpCuC^C^ measured by angle- 
resolved photoemission.EJ In order to reproduce the mea- 
sured spectral functions, one has to add longer-range hole 
hopping terms (f , t" and a-thce&,site hopping term) to 
the t-J model HamiltonianOO'E^I 

Motivated by the above discussion, we study the low 
energy states of the two-hole t-J model with different 
symmetries using exact diagonalization (ED). Previous 
ED studies were mostly performed on square lattices of 
16, 20, and 26 sites. Since we are focusing on the sym- 
metry of the wavefunctions, the geometry of the lattice 
is very important. For example, it is well-known that on 
the 16-site (4 x 4) lattice the system has spurious degen- 
eracy at the wavevectors (tt/2,tt/2) and (7r, 0). The 32- 
site lattice is free from these spurious effects. The d x 2_ y 2 
state of the two-hole model on the 32-site lattice has been 
studied in Ref. [?]. Besides, previous works have indicated 
that there is a p-wave state whose energy is very close 
to that of the d x 2_ y 2 state at realistic values of J/t. Ell 
The existence of this low energy state casts doubt on the 
robustness of the d x 2_ y 2 state as the ground state of the 
two-hole model. The energy difference is so small that 
any additional terms in the Hamiltonian may change the 
ground state from one to another. This cross-over may be 
realized by adding longer range hopping terms mention 
above to the t-J model Hamiltonian. Furthermore,|it has 
been shown that such terms favor the p-wave state.el This 
shows that the physics of the low-lying p state may not 
be irrelevant. Therefore it is important to understand 
the similarities and differences between the properties p£ 
the d x 2_ y 2 and the p states. In a previous publicationtll 
we have shown that the hole and current correlations of 
these two states are very different. In this paper, we will 
discuss in more detail the low-lying states of the two-hole 
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model in a larger range of J/t. The properties of these 
states will be compared at the realistic value of J/t = 0.3. 
Finally, we will demonstrate the effects of longer-range 
hopping and short-range Coulomb repulsion on the sym- 
metry of the ground state. 



II. LOW-LYING STATES WITH DIFFERENT 
SYMMETRIES 

The Hamiltonian of the t-J model is 

n = -t (4cj* + H.c.) + J^(Si ■ S, - ^mnj), 

(1) 

where <? and c are the projected fermion operators, and 
rii = cjcj is the fermion number operator. We solve this 
model with two holes on a 32-site lattice using the stan- 
dard Lanczos algorithm. To find the low-energy states, 
we consider subspaces with different rotational symme- 
tries and momenta. We concentrate on the rotational 
symmetries s, p, and d. States having larger angular mo- 
mentum usually have higher energy and are not relevant 
in our study. The momenta of the subspaces we study 
are (0,0), (7r,7r) and (7r,0). The dimensions of the sub- 
spaces with s and d symmetries are about 150 million, 
and that of the subspaces with p symmetry are about 
300 million. Most of the calculations are performed on a 
cluster of commodity personal computers. In subspaces 
with d symmetry, we achieve a speed of about 4.3 min- 
utes per Lanczos iteration on a 32-node cluster, where 
each node utilizes a 1 GHz Althon CPU. 

The energies of the low-lying states are tabulated in 
Table §. We are interested in the two-hole binding energy 
which is defined as 

Eb = E 2 h — ZEih + Eoh, (2) 

where E n h is the energy of the n-hole system. The ground 
state of the undoped system with energy E$h is a to- 
tally symmetric state with momentum (0,0). The ground 
state of the one-hole system with energy E\h has momen- 
tum (7r/2, 7r/2). Eb is the energy gain of the two-hole 
system relative to two one-hole systems, and therefore 
indicates the relative tendency for the two holes to form 
a bound pair. Keeping in mind that the "realistic" value 
of J is roughly between 0.2r; and OAt, we calculate the 
binding energies in a larger range of J in order to get 
a clear picture of the crossing-over of the energy levels. 
Fig. shows Eb of the low- lying states with s, p and 
d symmetries. The lowest energy state in the subspace 
with s symmetry is a spin singlet with momentum (0, 0). 
We call this the s state. At small J (J < O.lt), this is the 
ground state. But its energy rises very fast with J and 
soon changes from the ground state to an excited state. 
This shows that the s channel is energetically unfavor- 
able and is irrelevant at the realistic parameter J ~ 0.3t. 
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FIG. 1: Two-hole binding energies of the s (solid diamonds), d 
(solid circles), p^,-*) (solid triangles), £>(•„-, o) (empty triangles) , 
and P(o.o) (empty diamonds) states as defined in the text. 



The lowest energy state with d symmetry is a spin sin- 
glet with momentum (0,0) and d x 2_ y 2 symmetry. We 
call this the d state. This is the same state reported 
in Ref. f|. At J > 0.3i, it is the ground state of the 
two hole system. Note that the binding energy of this 
state scales linearly with J when J > 0.3i, reflecting the 
magnetic nature of the hole binding mechanism. E.2 Those 
states with p symmetry are more complicated. In Fig. |l| 
we show two low energy states with p symmetry. Both 
of them are two-fold degenerate and are spin singlets. 
One has momentum (tt,tt). We call it the p^^) state. 
The other has momentum (7r,0). We call it the p^.o) 
state.li^ In addition, we also show a higher energy p state 
with zero momentum, which we call the P(o,o) state. This 
state was mentioned in a previous study using a smaller 
lattice.Ej Also note that it is a triplet state. Fig. |l] indi- 
cates a general trend that as J decreases, the energy of 
the d state is pushed up faster than the p states. While 
at J > 0.3i the pu^) and p^.o) states are clearly excited 
states, the energy levels look quite complicated in the 
range 0.2t < J < 0.3i. Note that the d, p^,-*) an d P(tt,o) 
states are almost degenerate at J = 0.3i, although the d 
state is still the ground state. When J is slightly smaller 
than 0.3t, it looks like the p^.n) state has lower energy 
than the d state but at J = 0.2t, the p^o) state has the 
lowest energy among the three. At smaller J, we antic- 
ipate more serious finite-size effects. Based on our data 
only we cannot conclude where the crossing-overs of the 
energy levels occur, nor whether they exist in the ther- 
modynamic limit. But the fact that the crossing-overs in 
our system occur at J values within the realistic range is 
sufficient to cast doubt on the symmetry of the ground 
state in the two-hole model. Regardless of whether one of 
the p states actually becomes the ground state at certain 
range of J, their existence as low-lying states is unques- 
tionable. Furthermore, since the energy levels of these 
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TABLE I: Energies {E2h/t) of the low- lying states of the two-hole t-J model on a 32-site lattice at different J ft. The s, d, 
P(o,o); P(tt,tt) and P(tt,o) states are defined in the text. Also shown are the ground state energies Eih/t of the one-hole model. 
The ground state energy of the undoped model is Eoh ~ —37.7657342 J. 

J/t E lh /t E 2h /t 







s 


d 


P(0,0) 


P(tt,tv) 


P(t,0) 


0.1 


-6.419382 


-9.047472 


-9.015703 


-9.0208259142 


-9.017720 


-9.042840 


0.2 


-9.766819 


-11.903649 


-11.960871 


-11.9412605733 


-11.976890 


-11.983860 


0.3 


-13.161933 


-14.840663 


-15.045603 


-14.9685556227 


-15.045602 


-15.036724 


0.4 


-16.583058 




-18.189385 


-18.0477851139 


-18.163935 


-18.146186 


0.5 


-20.021272 




-21.368517 


-21.1594455768 


-21.312775 


-21.292519 


0.6 


-23.471807 




-24.572086 


-24.2941123207 


-24.483247 


-24.465752 


0.7 


-26.931796 




-27.793984 




-27.670397 


-27.660030 


0.8 


-30.399375 




-31.030391 




-30.871142 


-30.871578 



states are very close, inclusion of farther-than-nearest- 
neighbor hopping terms or Coulomb repulsion may easily 
change the ground state from one to another. Therefore 
we think that they are relevant and warrant detail stud- 
ies. In the following sections we study in detail the prop- 
erties of these states at J = 0.3i where they are almost 
degenerate. 



III. REAL SPACE STRUCTURE OF THE HOLE 
PAIR 
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Despite the fact that the binding energies of the d and 
p states at J = 0.3t are very similar and negative, their 
hole-hole correlations are very different. In Fig. we plot 
the hole-hole correlation function C(r) = ((1 — no)(l — 
n r )). As already discussed in Ref. |7] and also shown in 
Fig. ||, the holes in the d state attract each other, with 
C(r) having a maximum at r = \J~2. This has been oh-. 
served in previous ED studies using smaller lattices BlI 
Note that the observed attraction cannot be considered 
to be strong. The probability P(r) of the holes being at 
V2 [P(V2) = 0.23476] and V5 [P{V5) = 0.25037] are 
almost the same.E-1 The root-mean-square separation of 
the hole pair is r rms = yj (r 2 ) = 2.05865. Compared to 
the r rms of two uncorrelated holes in this lattice, which is 
2.38273, the hole pair in the d state is not tightly bound. 
This is consistent with its barely negative E^. On the 
contrary the holes in both p states seem to be mutually 
repulsive despite the negative E^. Their r rms are 2.53755 
and 2.53500 respectively, which are significantly larger 
than that of the d state. This shows that the hole pairs 
in the p^^) and P(-k,q) states are unbound. 

The spatial correlation of the hole pairs described 
above can be understood in terms of the nearest neigh- 
bor spin correlation (Si • Sj) at fixed hole configuration. 
This analysis was first carried out on a two-hole bound 
state using J.he technique of density matrix renormaliza- 
tion groupja and later by ED.El Fig. || shows our result 
for the d state on a 32-site lattice. The most intriguing 
feature is the existence of a strong next-nearest neighbor 
(i.e. across the diagonal of a square plaquette) singlet 
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FIG. 2: Hole-hole correlation function C(r) of the d and p 
states at J = 0.3t. Due to the lack of the rotational symmetry 
d (i.e., a rotation by n/2) in the p states, there may be two 
inequivalent points at the same r. In this case their values are 
indicated by the symbols and the solid line joins their mean 
values. 



bond in between two holes when they are at distance \/2 
apart. This is shown in Fig. ||a. It has been argueda that 
such singlet bond increases the hopping overlap between 
different hole configurations - when a hole hops this sin- 
glet bond becomes a strong nearest neighbor bond, and 
therefore is energetically favorable. But such singlet bond 
between two spins in the same sublattice causes spin frus- 
tration around the plaquette. Therefore the two holes 
tend to stay at next-nearest-neighbor distance in order 
to share the frustrating singlet bond. This explains why 
the holes are more likely to be found at distance \[2 as 
shown in Fig. @. When the holes are farther part as 
shown in Fig. 0b, a diagonal singlet bond still exists in 
the immediate vicinity of each hole, although it is weak- 
ened compared to the one in Fig. 0a. Besides, across 
each hole there is another strong singlet bond between 
two spins at a distance of 2. These two singlet bonds 
encourage the holes to hop towards each other because 
doing so will create two strong nearest-neighbor singlet 
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bonds. Such characteristic is retained in the other hole 
configurations shown in Fig. ||, although the strength of 
the singlet bonds weaken quickly as the holes move far- 
ther apart. Note that this set of hole configurations con- 
sistently have larger hole correlation C(r) as shown in 
Fig. |H For the other hole configurations where the holes 
are along the x or y directions, the characteristics of the 
spin correlations around the holes are very different. Two 
of these hole configurations are shown in Fig. H. When 
the holes are at nearest neighbor distance (Fig. pi), there 
are two strong nearest neighbor singlet bonds immediate 
above and below the hole pair. When one hole hops to 
become the hole configuration in Fig. |^a, one of these 
singlet bonds becomes the strong diagonal single bonds, 
and the other remains as a strong nearest neighbor bond. 
We can see the reminiscence of this bond on the four 
sides of the diagonal bond in Fig. ||a. Such hopping is 
strongly preferred as shown in the kinetic energy plot of 
Fig. [5| When the holes are at \[2 apart, the hopping 
energy along the four directions are significantly larger. 
When the holes are far apart (Fig. |]b), the spin correla- 
tion around the hole does not seem to possess any special 
feature. The hole hopping energy also decreases rapidly. 
However, contrary to a previous studyfH when the holes 
are close to each other, we do not see evidence for dimer- 
ization in the vicinity of the holes. In Fig. [j]we show the 
spin correlation relative to the undoped value. Except 
at the immediate vicinity of the holes, the spin correla- 
tion is rather unaffected and there is no clear evidence of 
dimerization. 

Next let us consider the ow-energy state Pi^^y As 
already pointed out beforep the spin correlation in the 
vicinity of a pair of near-by holes discussed above pro- 
vides a mechanism for hole binding. It has also been ar- 
gued that such spin correlation, is a result of the d x 2_ y 2 
symmetry of the wavefunctionB Then it is reasonable to 
expect that in wavefunctions with other symmetries, such 
mechanism for hole binding is missing and the holes may 
be unbound. Fig. seems to support this. When we com- 
pare the d and wavefunctions in the same subspace 

where the holes are at \[2 apart (Fig. ||a and ^a), we find 
striking differences in the spin correlation. The strong 
diagonal singlet bond between the holes found in the d 
state does not exist in the P( Wj7r ) state. The binding mech- 
anism discussed above is therefore missing in the p^.-n) 
state. This causes the probability for the hole pairs at \J~2 
apart to be small, as evident from Fig. ||. On the other 
hand, when the holes are far apart we observe a strong 
singlet spin pair across each hole (Fig. f7]b). It becomes a 
strong nearest neighbor singlet bond when the hole hops, 
thus increasing the overlap with other hole configurations 
and making it a preferred hole configuration. This fea- 
ture is missing in the d state (Fig. ||d) , where large hole 
separations are not preferred. This explains the repulsive 
nature of the holes in the P(tt,tt) state as shown in Fig. |^. 
Note that this interpretation is also consistent with the 
hole hopping energy as shown in Fig. [5]b. The hopping 




FIG. 3: Spin correlation (Si • Sj) in the d state at J = 0.3t. 
The two empty lattice points in each diagram are the loca- 
tions of the holes. Only nearest neighbor spin correlations are 
shown except in the vicinity of the holes where the correla- 
tions between spins at v2 and 2 apart are also shown. The 
width of the line joining two spins i and j is proportional to 
{St ■ Sj}/P(r), where P(r) is the probability of finding the 
holes at the relative locations as shown. Solid and broken 
lines mean {Si ■ Sj) are negative and positive respectively. 
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FIG. 4: Same as Fig. H but for different hole configurations. 
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FIG. 5: The hole hopping energy {c}cj) when the other hole 
is fixed at the location indicated by an empty circle, in (a) 
the d state and (b) the P(^ j7r ) state, both at J = 0.3t. Only 
a quarter of the lattice is shown. Values at other nearest 
neighbor pairs (ij) can be found by using the symmetry of 
the lattice. 
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FIG. 6: The spin correlation relative to the undoped value, 
(Si • Sj)/P(r) — (Si ■ Si) mldo ocd, of the d state at J = 0.3t. 



(Si ■ Sj} un dop C d = —0.34009. Solid and broken lines mean that 
the singlet correlation are enhanced (more negative) and sup- 
pressed (less negative) respectively relative to (Si • Sj) un dopcd- 
Only nearest neighbor correlations are shown. Those thick 
shaded lines mean that their values are too small to be shown 
in this scale, with magnitude smaller than 0.02. Their thick- 
ness do not represent the magnitude. They are drawn to show 
the sign of their values only. 
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FIG. 7: Same as Fig. g but for the P(^ j7r ) state at J = 0.3t. 
The scale of the line width is also the same as in Fig. [| 



energy is larger when the holes are farther apart. 



IV. MOMENTUM SPACE DISTRIBUTION 



Fig. ^| shows the electron momentum distribution func- 
tion (EMDF), defined as (flkcr) = ('^ker^ka') ' thrCC 

low energy states. Since in the two-hole model (nkf) = 
(n^i), we will leave out the spin index a. This quantity 
demonstrates another striking, difference between the d 
and p states. Previous worksta have shown that the gen- 
eral shape of (rik) is irrelevant to the structure of the hole 
pair. This is evident from the fact that in the two-hole 
model, An = ((rimo)) ~ { n (-K.ir))) is roughly the same as 
(An\ h + An} /l )o Note that this is true for all the three 
low-energy states. Nevertheless, those k along the mag- 
netic Brillouin zone boundary [(7r, 0) to (0, n)] are impor- 
tant in reflecting the structure of the bound state.Q The 
deviation of the EMDF from the undoped value, (rik) — \ , 
represents the hole weight at that k. As already pointed 
out in Ref. and shown in Fig. ||, (rik) of the d state 
along the magnetic Brillouin zone boundary has a max- 
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FIG. 8: EMDF of the low energy states. Each box represents 
an allowed momentum k. The three numbers inside each box, 
from top to bottom, are (rik) of the d, p^^) and P(n,o) states 
at J = 0.3t. Only the first quadrant of the Brillouin zone is 
shown, (ilk) in other quadrants can be constructed by using 
the symmetry of the Brillouin zone. 



imum at (tt/2, tt/2), and minimum somewhere between 
(7r/2,7r/2) and (tt,0). Note that in the d state (rik) is 
symmetric about the line from (0, 0) to (ir, 7r), and in the 
32-site lattice the two minima in the first quadrant of the 
Brillouin zone are at (37r/4, 7r/4) and (ir/4, 3tt/4). On 
the contrary, in the p states the minimum of (rik) along 
the magnetic Brillouin zone boundary is at (7r/2,7r/2). 
Note that the symmetry of the d state requires the hole 
weight at (7r/2,7r/2) to be zero. Hence its (nr. K /2,iv/2)) 
deviates from the undoped value of | by a minimal value 
of 0.0064. Such restriction is lifted in the p states and 
in fact (rik) shows that their hole weight is a maximum 
at (7r/2,7r/2). This feature resembles a hole pocket at 
(tt/2, tt/2). The fact that the EMDF of the p states have 
dimples at momenta corresponding to the single-particle 
ground state leads us to think that some form of rigid 
band filling approximation works in the p states but not 
in the d state. This can be accounted for by their differ- 
ent hole-hole correlations in real space. When the holes 
are farther apart, the overlap in the spin distortions they 
produce is smaller. Therefore the state can be better ap- 
proximated by some combination of the single-hole quasi- 
particle states, thus making the rigid band filling model 
a better approximation. 



V. SPIN CORRELATION 

Next we turn to the spin properties of the three low 
energy states at J = 0.3t . We first study the static 
structure factor S'(k) and the spin correlation function 
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FIG. 9: 
5(k). 
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Same as Fig. N but for the static structure factor 



(So • S r ). They are related by 

s(k) = J2 



— ik-r 



[So • S r 



(3) 



Fig. ^ shows the static structure factors. All three states 
show strong characteristics of a Neel state - 5(k) strongly 
peaks at (tt, tt). It is interesting to note that while those 
properties associated with the holes are very different 
between the p and d states, their static structure factors 
are very similar. As far as the spin property is concerned, 
the effect of the holes added to the system is to weaken 
the Neel order only, with no qualitative change. This is 
also indicated in the spin correlation function (So • S r ) in 
Fig. [To]. It is obvious that the spin correlations of the d 
and P(tt,tt) states have very similar behavior. If we fit a 
power law to the data points at r > 2, we find that the 
spin correlation decays as r -0 61 in the d state and r~ - 68 
in the p^^) state. They are to be compared to the spin 
correlation of the undoped ground state on the same 32- 
site lattice which decays as r~ - 25 . The spin correlation 
of the P( Tj o) state seems to be of a bit longer range and 
cannot be fitted satisfactorily by a power law. However, 
the difference does not show up significantly in the static 
structure factor. 

In the semi-classical theory of Shraiman and Siggia, a 
mobile hole produces long-range dipolar distortion in the 
spin background]^ A consequence is that the holes tend 
to stay away to minimize the overlap of the spin distor- 
tion they produce, resulting in dimples in the EMDF at 
momenta corresponding to the single-hole ground state. 
As discussed in the last paragraph of the previous sec- 
tion, the p states seem to fit this scenario better than 
the d state. Another consequence of the semi-classical 
theory is that in the presence of a small number of holos- 
the system may be unstable towards the spiral phase.Eil 
Therefore it is interesting to compare the tendency of the 
d and p states to have spiral order. A consequence of the 
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FIG. 10: Spin correlation function (So • S r ) of the three low 
energy states at J = 0.3£. Also shown is the same result 
on the ground state of the undoped system on the same lat- 
tice. Filled and empty circles represent positive and negative 
correlations respectively. 



spiral phase is that the maximum of S(k) will shift away 
from k = (tt, tt). As shown in Fig. ^, we do not see such 
pattern in any of the three low-energy states. However, 
the shift in k may be too small to be detected in our finite 
lattice with discrete k. In previous ED studies of the t-J 
model on smaller lattices, the shift in S(k) rwasi not ob- 
served until the doping level is much larger. ejI2j On the 
other hand, experiments has demonstrated incommen- 
surability in the magnetic fluctuation.cil Therefore one 
should not look for long range order that exists in the 
thermodynamic limit, but should instead look for short- 
range spiral order. For this reason we calculate the twist 
order parameter which is defined as 



X \ = S, x (S 4+x + S 



i+y) 



(4) 



Fig. [Tl] shows (xo ■ Xr) m the three low energy states. 
The short range behavior of this correlation function is 
very similar in the d and P( Vj7r ) states. In the p^.o) 
state it seems to decay faster. We find no enhance- 
ment in the short range spiral order in the p states rel- 
ative to the d state, and the correlation functions are 
very small in all three states. The susceptibility, defined 
as <7^IErX*| 2 ), are 0.02519, 0.02296, and 0.02492 for 
the d, p^.n) and p^o) states respectively. These are to 
be compared to the corresponding value of the undoped 
model on the same lattice, 0.02052. 



VI. ROBUSTNESS OF THE d x 
STATE 



,,2 GROUND 



Fig. |T] shows that in our system the energies of the 
d and p states are very close at J ~ 0.3t. This casts 
doubt on whether the d state is always the ground state 
of the two-hole model in the physically relevant range of 
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FIG. 11: Same as Fig. O except for the spiral spin correlation 
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J. Furthermore, the small energy difference means that 
including additional terms in the t-J model Hamiltonian 
may have profound effect on the selection of the ground 
state. This raises the question on the robustness of the d 
state as the ground state. In this section we discuss the 
effect on the energy levels when two kinds of terms are 
added to the t-J model. 

As pointed out in Ref. ||, a small size bound pair may 
be destroyed by a realistic short-range Coulomb repul- 
sion. On a two-leg ladder, the authors have shown that 
the bound pair disappears around V ~ 4 J, where V is the 
nearest neighbor Coulomb repulsion. Recently Coulomb 
repulsions are found to enhance the staggered orbital cur- 
rent on a two-leg ladder.Ba On a two-dimensional lat- 
tice, such staggered curr.eM|U|has been found even with- 
out Coulomb repulsions S3c3 Adding them to the two- 
dimensional model may do more than just pushing the 
holes apart. The repulsions should weaken hole bind- 
ing and raise the energy of the d state. But in the p 
states, the effect is less significant because the holes are 
already far apart. As a result one would expect short- 
range Coulomb repulsions to favor the p states to the d 
state. Here we add two repulsion terms V\ and V% to 
the t-J Hamiltonian. They are repulsions between two 
holes at 1 and v2 apart respectively. To test the sta- 
bility of the d x 2_ y 2 ground state, we choose to use the 
values Vi — 0.6t, Vi — 0.3t. The binding energies of the 
three states are shown in Fig. [l^ as a function of J/t. As 
expected, the binding energy of the d state is pushed up 
more than that of the p states. At J = 0.3t, only Eh of 
the P(7r,7r) state is barely negative. In the range we study, 
Q.2t < J < 0.8, the p^.-k) state has the lowest energy. 
This calculation shows that even a small Coulomb repul- 
sion will destroy the two-hole bound state and select the 
p state as the ground state. 

It is know that in order to reproduce the single hole dis- 
persion observed experimentallyjla one has to add Ion err, 
range hole hopping terms t' and t" to the t-J modeliiSt^l 



TABLE II: Energies Eik and binding energies Eb of the lowest 
energy states with different symmetries of the two-hole t-t'- 
t"-J model with J = 0.3t, t' = -0M, and t" = 0.2t Ground 
state energy E^h of the one- hole model is — 13.716526t. 



symmetry 


momentum 


spin 


E 2h /t 


E b /t 


d 


(0, 0) 





-15.588993 


0.51434 


s 


(0, 0) 





-16.074058 


0.02927 


V 


(0, 0) 


1 


-16.080823 


0.02251 


V 


(tt,0), (0,tt) 





-16.109653 


-0.00632 


p 


(VT, 7t) 





-16.164444 


-0.06111 



These terms are relevant here for two reasons. First the. 
magnitude of these terms, t' — — 0.3t and t" — 0.2tf3 
are larger than the energy difference between the d and 
p states at J = 0.3t. Therefore they may have serious 
effect on the low- lying energy levels. Second, it has been 
argued that the form of single-hole dispersion resulting 
from the t' and t" terms favors the p statefl In Ref. [l^, 
three-site hopping terms are also included in addition to 
t' and t" in order to reproduce the spectral functions 
measured experimentally. However, if we just want to 
reproduce the single hole dispersion, it will be suffice to 
include the t' and t" terms only. Table [II] shows the low- 
est energy state in different symmetry subspaces when 
t' = — 0.3i and t" = 0.2t are included in the t-J model at 
J = 0.3t. We can see that the structure of the low-lying 
energy levels is seriously affected. First, the energy of the 
d state is pushed up by a lot compared to other states. 
Its energy is now even higher than an s and a (triplet) p 
state. Note that in the t-J model these two states have 
higher energies than the d state (see section II). Nev- 
ertheless, their binding energies are positive. Only the 
two singlet p states have negative binding energies. The 
p state with momentum (tt, if) now becomes the ground 
state. Note that its binding energy does not differ too 
much from that of the t-J model as shown in Fig. [I], 
which is -0.05146*. 



VII. CONCLUSION 

We tried to map out the low-lying states of the t-J 
model with two holes on a 32-site square lattice at vari- 
ous J/t. While previous studies have focussed on finding 
the critical J/t where hole binding is lost, we concen- 
trate on the symmetry of the ground state at different 
J/t. We find that at large J/t, the ground state has 
d x 2_ y 2 symmetry as reported before. But there are low- 
lying states with p symmetry whose energies come closer 
to the ground state as J/t decreases. At J/t < 0.3 one 
of these p-wave state becomes the ground state. We note 
that these p-wave states have non-zero momenta. Pre- 
vious numerical works have revealed similar finite mo- 
mentum states. But calculations carried out on lattices 
smaller than 32 sites are biased by their geometry. As 
a result, these finite momentum p states have not re- 
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FIG. 12: Same as Fig. [l] except that the Coulomb repulsion 
terms V\ = 0.6t and V2 = 0.3t are added to the t-J model 
Hamiltonian. The states shown are d (circles), P(^ l7r ) (solid 
triangles) and P(^-.o) (empty triangles). 

ceived much attention. Nevertheless, a similar cross-over 
in the symmetry of the ground state from d to p sym- 
metry has been, reported in a series expansion study of 
the t-J modelrH There the cross-over was found to be at 
J/t ~ 0.4, which is quite close to our value of J/t ~ 0.3. 
This crospjieyer was also observed in the anisotropic t- 
J z model.ul23 Although finite-size effects prevent us from 
concluding definitely whether this cross-over exists in the 
thermodynamic limit, the fact that it occurs within the 
parameter range of interest make us feel that the p states 
are relevant to the low energy physics of the t-J model. 
Further investigation reveals qualitative differences and 
similarities in the properties of the p and d states. Per- 
haps the most intriguing is the difference in their hole 
properties. The holes are mutually attractive in the d 



state, but are repulsive in the p states. The spin struc- 
ture in the vicinity of the holes provides an intuitive ex- 
planation of the contrasting behaviors. And as a result of 
this, the electron momentum distribution functions are 
also very different. While the EMDF of the d state shows 
no sign of hole pockets, that of the p states clearly have 
dimples at (tt/2,tt/2). Such hole properties suggest that 
the rigid-band filling assumption should work better in 
the p than in the d state, if it works at_all. And in the 
framework of the semi-classical theory]^ 1 ] one would ex- 
pect the p states to have stronger tendency to show spiral 
spin order than the d state. However, we do not find any 
enhancement in the short-range spiral spin correlation in 
the p states over the d state. In fact the spin correlations 
of the d and p states are very similar and show signals of 
the Neel state. 

We have also demonstrated that the d x 2_ y 2 ground 
state is not robust. It can be destroyed easily by includ- 
ing realistic terms to the t-J model Hamiltonian. Both 
the short-range Coulomb repulsion and longer range hop- 
ping terms (t' and t") favor the p state with momentum 
(-7T, 7r) as the ground state. This shows that hole pairing 
in the d x i_ y i channel may not be a generic feature of the 
t-J model. There is a competing p state which has no 
hole pairing. The implication of our result is that the 
symmetry of the two- hole ground state of the t-J model 
is yet to be determined. This is in contrary to a recent 
numerical studyEj which concludes that pairing in the 
d x 2_ y 2 channel is a robust property of the t-J model. 
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